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We develop an analytical technique to derive explicit forms of thermodynamical quantities within 
the asymptotic approach to non-extensive quantum distribution functions. Using it, we find an 
expression for the number of particles in a boson system which we compare with other approximate 
scheme (i.e. factorization approach), and with the recently obtained exact result. To do this, we 
investigate the predictions on Bose-Einstein condensation and the blackbody radiation. We find that 
both approximation techniques give results similar to (up to 0[q — 1)) the exact ones, making them 
a useful tool for computations. Because of the simplicity of the factorization approach formulae, 
it appears that this is the easiest way to handle with physical systems which might exhibit slight 
deviations from extensivity. 
PACS Number(s): 05.20.-y, 05.30. Jp, 05.30.Fk 

I. INTRODUCTION 

Since the papers by Tsallis non-extensive statistical formalism has been shown to be not only robust -ii allows 
generalizations of all necessary fundamental concepts of thermostatistics but also useful -it provides a suitable 
theoretical tool to explain some of the experimental situations where standard thermostatistics has shortcomings, due 
to the presence of long-range interactions, or long-range memory effects, or (multi)-fractal space-time constraints-. 
See Ref. Q for a periodically updated bibliography. 

The core of this generalized formalism is defined through a generalized entropy 

= k ^_^='^\ {q e 5R), (1) 

where A; is a positive constant, {pi} is a set of probabilities and W is the total number of microscopic configurations. It 
is easy to verify that the q ^ 1 limit immediately recovers the usual (extensive) Boltzmann-Gibbs entropy. Moreover, 
if a composed system A + B has probabilities which factorize into those corresponding to the subsystems A and B, 
then Sq{A + B)/k = Sq{A)/k + Sq{B)/k + (1 - q)Sq{A)Sq{B)/k'^. This property clearly exhibits the fact that the 
parameter q characterizes the degree of non-extensivity of any physical system. 

The generalization of quantum statistics for non-extensive systems was only accomplished, up to recent days, in an 

a approximate fashion, by using two different schemes. One of them is the Asymptotic Approach (AA), of Tsallis et al. 
the other one is the Factorization Approach (FA), of Biiyiikkilig et al. The physical applications studied so 
far within these two approximations include the blackbody radiation the Stefan-Boltzmann constant P-pd]|, 

C5 ' and some aspects of the early universe physics Moreover, the AA has also been used in some other works 

O such as the Bose-Einstein condensation ||l^ , the specific heat of ^He [|5| , thermalization of an electron-phonon system 
@ and cosmology 0,|l8|. Although some detailed analysis on these approximate schemes suggest that both 
schemes could be helpful in physical applications -at least, for (1 — g)-order corrections-, this was still doubtful. A 
complete verification needed a comparison between the results of these approximate schemes and the exact ones. But 
an exact treatment of non-extensive quantum distributions was not available up to the recent papers of Rajagopal 
\ et al. [|2^,^. Just after this work, Lenzi and Mendes have also given an exact treatment of blackbody radiation 
p2[ . All these recent efforts enable us to make a comparison between the approximate and exact schemes, which will 
ultimately show whether the AA and the FA are useful or not. This will be the main purpose of this paper. 
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In Section II, we review the approximate and exact results and develop an analytical method to derive the explicit 
form of any measurable quantity within the AA. We compare the approximate and exact results in Section III using 
(i) the predictions of one of the experimental tests suggested in Q and (ii) the blackbody radiation. Finally, we give 
our final comments in Section IV. 



II. NON-EXTENSIVE QUANTUM STATISTICS 
A. Asymptotic approach 

Within the AA, namely in the /3(1 — g) — > 1 limit, the generalized partition function is given by 

Z,c^zJl-l{l-q)/3^ (h 



from where the generalized distribution function of non-interacting bosons can be found, up to (1 — (7)-order, as 



(n), = (n), + {l-q) (n), { ln(l/Zi) + [x - ^Jj) 



(«^)r~ 



where x = (3e , ip = Pfi, and 
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aX — lp 



The standard (g = 1) partition function is given by 

Zi = 



1 — e-i^-^) ' 



(2) 



(3) 



(4) 



(5) 



This approximation has found a wide range of applications up to now, however, no attempt has been made for deriving 
some of the thermodynamical quantities within this approach, directly using Eq. (^). 

One aim of this paper is to provide a technique for computing, in a closed form, the kind of integrals needed to find 
the average number of particles within the AA. To do this, let us start by writing down the definition of the average 
number of particles: 



27ry(2mfc)3/2r3/2 



(6) 



where all variables have the usual meaning. Using Eq. (^) and the definitions of x and ip^ this expression turns out 
to be 



'—j^ [1st + (1 - g)(/2 + h)] , 



where 



and 



x^^'^dx 



(x — ^)x^^'^dx 



h = 



{x — ip)x^/'^dx 
e^-'^ - 1 



(7) 



(8) 



(9) 



I2 and /s are the {q — 1) order correction to the standard {q — 1) result and here 1st stands for the standard integral 
appearing in the solution of the extensive case pSl]. 1st and I2 have standard forms, and could easily be solved as: 
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Ist=T{i/2)g^/2{z), 



(10) 



2 



oo 5/2-1 J j-oo 3/2-1 J 

^j^^ = r(5/2).g3/2(^)-^r(3/2)g3/2(^), (h) 



^ "'0 



where z is the fugacity and is defined as z = e*^^. On the other hand, is more involved, and it takes the form: 

h = a + b-ic- d, (12) 

where 

^= / r ,.2 ' (13) 



[e==-'/' - 1] Jo [e^-"^ - 1] 

(a; - iifx^l'^e'-^dx , /■°° (x - V)2a;i/2e2(^-V')da; 



[e^-'/' - 1]^ ' io [e^-'^ - 1]^ 



(14) 



In an Appendix, we provide an analytical technique (maybe there are others) to compute each one of these integrals. 
Using this technique, we obtain the average number of particles as: 

2.F(2l'fcr)3/2 - r(3/2),3/2(.) + (. - 1)V^ X 

"?! Q 7 1 . q . 1 

(15) 



3 9 7 1 9 

293/2(2:) - -35/2(2) + -^^53/2(2) - 2-091/2(2:) - 2V'^5l/2(2:) + ^-0^5-1/2(2:) 



Here, (J^n) ^ stands for the number of particles in the excited states (e 7^ 0). As in the standard case, we have separated 
the contribution of the state given by e 0, which has zero weight in the integrals. For this level of energy, we found. 



W,(£ = 0)^^^l + (g-l) 



zlnz Inz z(lnz)^ (Inz) 
lnz + +3- I ' ^ > 



1-z 1-z (1-2)2 [\-zf 



(16) 



When z ^ 1, all the correction terms go to zero. When z ^ 1, some of the terms are divergent but the usual shape 
is unchanged. This can be seen in Fig. 1. 

Numerical analysis, which we show in Fig. 2, illustrates that the maximum correction is attained for z = 1. Then, 
the number of particles in all excited states is bounded by, 

TO, <^^M^ [2.315 + (,-1,4271^ (17) 

It is worth noticing that the AA is such that not all terms in the (1 — q) correction are positive (or negative, depending 
on the choice of q) definite. Moreover, their maximum values are not always attained at middle points of the interval 
of interest, and although they do have bounded expressions, the maximum correction is obtained only for z = 1. This 
differs from what happened in the FA, where each term had a maximum value within the interval of interest [^ . 
The order of magnitude of the maximum correction is, however, the same in both approximations. 

Any interested reader could easily apply the same technique, which we introduced in the Appendix, to compute 
any other therniodynamical quantity, whenever it is needed. 

B. Factorization approach 

Within the FA ||^, the generahzed distribution function of bosons is given, at (1 — q) order, by Q 

(•7. _ ib'S^e^-'^ 

W,-(»)i + (g-l) ,, ^ ,,2 : (18) 
9 2(e:^-'''-l) 
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where (n)-^, x and 4' have the same definitions as before. At this point, the remarkably simpler form of this result, 
when compared to the result of the AA [Eq. (^)], is worth emphasizing. 

In this context, we have found general expressions for some thermodynamical quantities of bosons and fermions p4[ ; 
here we quote only the average number of particles for bosons, since it will be adequate for our proposed comparisonfH: 

(iV.),<?M^ [2.315 + ,,-1)3.0791. (19) 



C. The exact result 

Although the results of the A A and the FA have been successfully used in a wide range of physical applications, 
an exact treatment of non-extensive quantum statistics was lacking until the recent work of Rajagopal, Mendes and 
Lenzi pc|j2l[| . In their analysis, they have given the many-particle q-Green function in terms of a parametric contour 
integral over a kernel, multiplied by the usual grand canonical one particle Green function which now depends on q. 
They managed to obtain exact expressions for thermodynamical quantities, such as (N)^. 

To proceed further, let us quote here some of the results of pO| , ^ . Rajagopal et al. have used the general contour 
integral of the form, 

/ duexp{-bu){-u}-^ = -1^, (20) 
27r Jc r(z) 

with 6 > and Re z > 0, and where the contour C starts from -l-oo on the real axis, encircles the origin once 
counterclockwise and returns to -l-oo. Using the g-Green functions, and after some cumbersome algebra, they finally 
obtain (for bosons) 



Ic " J-oo 27r 7 (27r)^ [e-/3(i-9)"('^-A') - 1" 
where D is the dimension of space, A['p]uj) is the spectral weight function and 



A{p;lo), (21) 



KfHu) = *^^^^^|^exp(-.)(-.)-V(i-.), (22) 

and 

Z,(/3,Ai)= / duK'^^\u)Zi{-Pu{l-q),y). (23) 
Jc 

This exact expression for the average number of particles finally gives us the opportunity to make a comparison 
between the exact and the approximate results. 

III. EXACT AND APPROXIMATE RESULTS 

A. Bose-Einstein condensate 

One of possible experimental tests of the validity of the q-framework is based on a recent work on Bose-Einstein 
condensation of a small number of atoms (of the order of 100 to 170), confined to a small region of space by magnetic 
trapping |^^. By taking free particle spectral weight function, namely A{p;uj) — 2tt6{uj — fF/2m), near the Bose- 
Einstein condensation, they have found 



^We take advantage here to signal out a mistake in the last equation of Ref. where the correction appears to be proportional 
to 0.886 (q — 1) and should have a minus sign in front of it pa]. 
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{Tc, 



3/2 



{N)i C(5/2) 



(i-'z)^/2r(f3f + i 



C(3/2) 



3/2 



rffEf + 2 



p 2^ 

1-9 



p 2^ 

1-9 



(24) 



Here, the ~ sign reflects that this equation is vahd near the Bose-Einstein condensation, which does not change the 
fact that it is an exact result. We may now expand this expression in powers of {q — 1). Up to first order. 









3/2 1 



l + (g-l) 0.456- 0.023 (iV)i 



(Teh 



3/2N 



The two previous equations deviate from each other very soon when {q — 1)^ is not negligible. 

Let us now derive similar expressions for the AA and FA in order to compare them with Eqs. 
condition for the appearance of Bose-Einstein condensation can be expressed as 



(25) 



and (E5h. The 



(26) 



Alternatively, with constant (N)^ and V, using Eq. (17) for the AA and Eq. (19) for the FA, this condition can be 
recast in the form 



T < {T, 



2/3 



(27r)3/22mfc [ V [2.315 + (q - 1)k] 
where k — 4.27 for the AA and k = 3.078 for the FA. Then, we can organize these expressions to give 



3/2 



[2.315+ (q- 1)k] 
2.315 ■ 



(27) 



(28) 



Eq.(17) and (19) have corrections which are trivial (not depending on z) just because we have approximated them: 
the actual complete results are Eqs. (15) and (16) for the AA, while those for the FA can be found in our previous 
paper [24] . We managed the dependence on z in order to obtain an upper bound for the corrections and simplify the 
analysis that follows. Differences between Eqs. ( p8[ ) and ( p5| ) are worth noticing: the later depends on {Tc)q and (N)-^ 
in a much stronger way. However, as we shall see, for q close to 1 these differences are not important. 

We would now like to choose physically suitable q values. An early Universe test based on the FA has shown |l^] 
to produce a bound [g — 1| < 4.01 x 10^^^ thus we have q = 0.996. The other q value which we use comes from 
a very recent work on pion transverse-momentum correlations in Pb-Pb high-energy nuclear collisions [^ . In that 
work, a deviation of [g — 1| = 0.015 from the standard statistics is found to be sufficient for eliminating the puzzling 
discrepancy between theoretical calculations and experimental data |2^. Thus, we shall use q = 0.985 (in fact, in 
p7| , q = 1.015 has been used, but since the exact result is given for g < 1 values, we must take q = 0.985, which has 
the same [<? — 1| deviation). In Fig. 3 we plot (N)^ j {N)^ versus (Tc), / {Tc)i for two representative values of {N)^, 
and the two quoted values of q. However, note again that in our approximated schemes, (N)^/ (N)-^ as a function 
of {Tc), I (Tc)]^ is in fact independent of the particular value of (N)-^. From Fig. 3, the following conclusions can be 
drawn: At the order of such q values, the AA and the FA are almost the same, and in (1 — (7)-order correction, any 
of them could be used with the same confidence (maybe the FA would be preferable due to its remarkably simpler 
form). Only in those situations of extremely high experimental precision one could distinguish between the exact and 
approximate results. 



B. Blackbody radiation 

Very recently, an exact analysis of the blackbody radiation within the q- framework has been given [ p2[ . This exact 
analysis gives the generalization of the Stefan-Boltzmann law as 



5 



where 



.J ^ 5kT^ ^ m - g)/(l - g)] 

Z| ml r[(2 - q)/{l -q)+ 3(m + 1)] ' ^ 



r[(2-g)/(l-g)] 

" !r[(2-g)/(l-g) + 3m]' ^-""^ 



m=0 

and 



^ _ 4r(3)C(4) /27rV^/^kT 

~ [27ri/2(i_g)]3r(3/2) 
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Let us now recall the Stefan-Boltzmann law derived by using the AA [pUlCj] and the FA 



where 6 = 40.018 for the AA and = 62.215 for the FA. 



(31) 



U,^ ^3^3 [6.4939 -(l-g)g] (32) 



For the comparison of the exact and approximate Stefan-Boltzmann laws, we again try to choose a value of q which 
is in accordance with the blackbody radiation. The possible (/-correction could be at the order of 10"''' or 10~^. Thus, 
here we shall use again the largest deviation predicted for the g-correction namely I? — 1| < 5.3 x 10^^, which 
gives q — 0.99947 . In Fig. 4 we present the behaviour of the exact [Eq. (|2^)] and the approximate results [Eq. (|3^)] 
for q — 0.99947. It is seen from the figure that for such order of g-correction the approximate results are very close to 
the standard {q = 1) case without exhibiting any curvature, contrary to the exact result. 



IV. FINAL REMARKS 



We have managed to develop an analytical technique to express thermodynamical quantities for the asymptotic 
approach of quantum distribution functions. We have shown that, for simple boson systems, and for all g-values 
admitted by the existing bounds, both approximate schemes (the AA and the FA) are in agreement with the exact 
result (see figures) . The magnitude of the deviation is quantified in previous formulae and could be seen if there is 
enough experimental precision. Otherwise, the simpler form that the factorization approach exhibits makes a case for 
its use as a standard and safe procedure for (1 — g)-order corrections. 
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FIGURE CAPTIONS 



Figure 1 : Nq{e = 0) of the AA as a function of z for different q values. 

Figure 2 : (a) The contribution of the different terms that enter the {q — 1) correction to the average number of 
bosons within the AA. On the right corner of the figure, the curves corresponds to the following order: first, fourth, 
sixth, fifth, second and third term, (b) The total (q — 1) correction to the average number of bosons within the AA. 
Its maximum possible value is attained at z = 1, and the correction goes as 77^/^(3/2^(3/2) — 9/8^(5/2)) = 4.27 . 

Figure 3 : Bose-Einstein condensation: Plot of {N)^ / (N)-^ as a function of (Tc)^ / (T'c)j for (a) q = 0.996 and (b) 
q = 0.985 . 

Figure 4 : Blackbody radiation: Internal energy versus {2-jTkT) / (he) for q = 0.99947 . 
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APPENDIX A: 



The calculation of the integral a can be done as follows. Let us define, introducing an extra parameter m, 

Jo 

Then, we could write. 



a = 





'=Jo 






' dl' 




" {x-ijj)x^/^dx 




dm 


m— 1 


Jo {e^~^ — mf' 


m=l 



Defining m = e"^, it is easy to write, 



d 



I' 

Jo 



{x — tj))x^/'^dx 



dm Jo m{m~'^e^~^ — 1) dm 
where = tjj — (p. This let us to obtain, 

■ 1 



gX-V _ 1 



°° x^/^dx 



_d_ 
dm 



m 



r(5/2)55/2(V') - ^r(3/2)ff3/2(^')) 



Finally this gives us the solution of the integral a: 

a = r(5/2)ff3/2(^) - r(5/2)55/2(^) - Vr(3/2)5i/2(^) + ^r(3/2)ff3/2(^). 
For the calculation of c, we may conveniently define ai as, 

(x — ip)x^^'^dx 



ai 



and derive with respect to I to obtain 

dai 
'dl 



;=i 



[ei(x-V') _ 1]^ 
{x - i)fx^l'^e^-^dx 



= -2c. 



Now, to compute a;, we may change variables as follows: 
We then obtain, 



ai 



Deriving with respect to Z, we get 

5 



dai 

la ~ 2^72 



{r(5/2) [gs/2{^) - 55/2(V5)] + V;r(3/2) [g^,S) ' 51/2^] } + 

^ {r(5/2) [^5i/2(V^) - i^gmW] + ^r(3/2) [53/2^ - 5i/2(^)] 

+ ^r(3/2) [V'ffi/zl^) - V'ff-i/2(^)] } 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 



(A8) 



(A9) 



(AlO) 



From this equation, the integral c can be obtained by making 1 = 1 and i/j = ip. Here, we should note that, in all 
previous calculations, we have used (i) the result 



5„-i(.) = z^b„(z)] = ^M^)]. 



(All) 
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and (ii) Robinson's power series representation pq| (which is valid for all n values) : 



g„(a) = r(l - + ^ i-^C(ri - 0« 



(=0 



(A12) 



for the gn functions, where a = — ln(z) and Q is the Riemann zeta function. Using this result, one can recover the 
relationship for the derivatives of gn functions, and use it to evaluate g-i/2- We finally obtain the result for the 
integral c: 

c = ^r(5/2)53/2(^) - ^r(5/2).95/2 W + ^r(3/2)Vg3/2(V') - \n?,/2)i^g,/M) 
-ir(5/2)^5i/2W + ir(5/2)V'53/2(V') - ir(3/2)Vg3/2(V') + ir(3/2)V;ffi/2(V') 



-ir(3/2)V2^i/2W + ir(3/2)V'Vi/2W. 



(A13) 



We proceed further to compute integral h. To do so, we apply the following procedure. Let us define a new integral 
with an extra parameter, such that 



1 



,3/2 



1 



rr(3/2).g3/2(^), 



/o e^*"''') - 1 

where we have used the change of variable x = mx and ip = imp. If we now derive with respect to m, we obtain. 



(AM) 



dam 






dm 


m— 1 


. Jo 



(^gm(a:—i/') _ 1)2 



Since we have. 



,r(3/2)53/2(^) 



dm 2m^/^ 
it is easy to write down the solution of the integral b 



,3/2 



r(3/2)V5i/2(^), 



(A15) 



(A16) 



b = -r(3/2)53/2(^) - r(3/2)V'<7i/2(V')- 



(A17) 



Once the integral b is calculated, the integral d can be obtained as follows. Let us define Um as above. Deriving it 
twice with respect to the parameter m, we obtain 



d^a„ 



dm^ 



2d~ I„ 



where, 



^new 



(a; - -0)^a;i/2e'"("=~^)dx 



|^gm(a;-V) _ l] ^ 

It is easy to compute this integral with a similar trick. We need to define, with usual notation, 

(a; — 'il;)x^/^dx 1 



r(5/2)g5/2(^) - Vr(3/2)53/2(V') 



/o [e'(^-'A) - 1] ;5/2 
and derive it with respect to I. Further evaluation in Z = 1 reproduces Inew'- 

Ine^ = ^r(5/2)g5/2(V') " v4r(3/2).93/2(V') - V'r(5/2)(73/2(V') + V'r(3/2).9i/2(V'). 
Thus, we finally have the solution of the integral d: 

d = ^r(3/2)53/2(V') - ^mm9i/2W + Jv'r(3/2)5_i/2(^) + ^r(5/2)g5/2(V) - 
V^r(3/2)g3/2(V') - ^V^r(5/2)g3/2(^) + iv^r(3/2)53/2(^) + iv'r(3/2)5i/2(^). 



(A18) 



(A19) 



(A20) 



(A21) 



(A22) 



10 



// 



each term of ^-correction 




Fig 2(b) 




Fig 3(a) 




Fig 3(b) 
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